Search for: All records

Marine δ¹⁸O data reveal astronomical forcing of the climate and cryosphere during the Miocene, when atmosphericPco₂was on par with emissions scenarios over the next century. This inspired hypotheses for how Milankovitch cycles, ice-ocean interactions, and greenhouse gases influence ice volume. Mass balance controls for marine and terrestrial ice sheets differ, and proxy data collected far from Antarctica provide valuable but limited insight into regional processes. We evaluate clast abundance data from Antarctic marine sedimentary records, observing a strong signal of eccentricity and precession coincident with a terrestrial ice sheet and a clear obliquity signal at the margins of a marine ice sheet. These analyses are integrated with a synthesis of proxy data, and we argue that high variance in obliquity forcing (mediated and enhanced by the ocean and atmosphere) can inhibit ice sheet growth, even when insolation forcing is conducive to glaciation. This “obliquity disruption” explains cryosphere variability before the existence of large northern hemisphere ice sheets. 

Climate and ecosystem dynamics vary across timescales, but research into climate-driven vegetation dynamics usually focuses on singular timescales. We developed a spectral analysis–based approach that provides detailed estimates of the timescales at which vegetation tracks climate change, from 10¹to 10⁵years. We report dynamic similarity of vegetation and climate even at centennial frequencies (149⁻¹to 18,012⁻¹year⁻¹, that is, one cycle per 149 to 18,012 years). A breakpoint in vegetation turnover (797⁻¹year⁻¹) matches a breakpoint between stochastic and autocorrelated climate processes, suggesting that ecological dynamics are governed by climate across these frequencies. Heightened vegetation turnover at millennial frequencies (4650⁻¹year⁻¹) highlights the risk of abrupt responses to climate change, whereas vegetation-climate decoupling at frequencies >149⁻¹year⁻¹may indicate long-lasting consequences of anthropogenic climate change for ecosystem function and biodiversity. 

Abstract. Relating stratigraphic position to numerical time using age–depth models plays an important role in determining the rate and timing of geologic and environmental change throughout Earth history. Astrochronology uses the geologic record of astronomically derived oscillations in the rock record to measure the passage of time and has proven to be a valuable technique for developing age–depth models with high stratigraphic and temporal resolution. However, in the absence of anchoring dates, many astrochronologies float in numerical time. Anchoring these chronologies relies on radioisotope geochronology (e.g., U–Pb, 40Ar/39Ar), which produces high-precision (<±1 %), stratigraphically distributed point estimates of age. In this study, we present a new R package, astroBayes, for a Bayesian inversion of astrochronology and radioisotopic geochronology to derive age–depth models. Integrating both data types allows reduction in uncertainties related to interpolation between dated horizons and the resolution of subtle changes in sedimentation rate, especially when compared to existing Bayesian models that use a stochastic random walk to approximate sedimentation variability. The astroBayes inversion also incorporates prior information about sedimentation rate, superposition, and the presence or absence of major hiatuses. The resulting age–depth models preserve both the spatial resolution of floating astrochronologies and the accuracy as well as precision of modern radioisotopic geochronology. We test the astroBayes method using two synthetic datasets designed to mimic real-world stratigraphic sections. Model uncertainties are predominantly controlled by the precision of the radioisotopic dates and are relatively constant with depth while being significantly reduced relative to “dates-only” random walk models. Since the resulting age–depth models leverage both astrochronology and radioisotopic geochronology in a single statistical framework they can resolve ambiguities between the two chronometers. Finally, we present a case study of the Bridge Creek Limestone Member of the Greenhorn Formation where we refine the age of the Cenomanian–Turonian boundary, showing the strength of this approach when applied to deep-time chronostratigraphic questions. 

Paleohydrologic proxy data and climate models show how and why eccentricity and precession influenced early Eocene hydroclimate. 

Lacustrine strata are often among the highest-resolution terrestrial paleoclimate archives available. The manner in which climate signals are registered into lacustrine deposits varies, however, as a function of complex sedimentologic and diagenetic processes. The retrieval of reliable records of climatic forcing therefore requires a means of evaluating the potential influence of changing sedimentary transfer functions. Here, we use high-resolution X-ray fluorescence core scanning of the Wilkins Peak Member of the Green River Formation to characterize the long-term evolution of transfer functions in an ancient lacustrine record. Our analysis identifies a shift in the frequency distribution of Milankovitch-band variance between the lower and middle Wilkins Peak Member across a range of temporally calibrated elemental intensity records. Spectral analysis of the lower Wilkins Peak Member shows strong short eccentricity, obliquity, precession, and sub-Milankovitch−scale variability, while the middle Wilkins Peak Member shows strong eccentricity variability and reduced power at higher frequencies. This transition coincides with a dramatic decline in the number and volume of evaporite beds. We attribute this shift to a change in the Wilkins Peak Member depositional transfer function caused by evolving basin morphology, which directly influenced the preservation of bedded evaporite as the paleolake developed from a deeper, meromictic lake to a shallower, holomictic lake. The loss of bedded evaporite, combined with secondary evaporite growth, results in reduced obliquity- and precession-band power and enhanced eccentricity-band power in the stratigraphic record. These results underscore the need for careful integration of basin and depositional system history with cyclostratigraphic interpretation of the dominant astronomical signals preserved in the stratigraphic archive. 

The Wilkins Peak Member (WPM) of the Green River Formation in Wyoming, USA, comprises alternating lacustrine and alluvial strata that preserve a record of terrestrial climate during the early Eocene climatic optimum. We use a Bayesian framework to develop age-depth models for three sites, based on new 40Ar/39Ar sanidine and 206Pb/238U zircon ages from seven tuffs. The new models provide two- to ten-fold increases in temporal resolution compared to previous radioisotopic age models, confirming eccentricity-scale pacing of WPM facies, and permitting their direct comparison to astronomical solutions. Starting at ca. 51 Ma, the median ages for basin-wide flooding surfaces atop six successive alluvial marker beds coincide with short eccentricity maxima in the astronomical solutions. These eccentricity maxima have been associated with hyperthermal events recorded in marine strata during the early Eocene. WPM strata older than ca. 51 Ma do not exhibit a clear relationship to the eccentricity solutions, but accumulated 31%−35% more rapidly, suggesting that the influence of astronomical forcing on sedimentation was modulated by basin tectonics. Additional high-precision radioisotopic ages are needed to reduce the uncertainty of the Bayesian model, but this approach shows promise for unambiguous evaluation of the phase relationship between alluvial marker beds and theoretical eccentricity solutions. 

ABSTRACT The Green River Formation preserves an extraordinary archive of terrestrial paleoclimate during the Early Eocene Climatic Optimum (EECO; ∼ 53–50 Ma), expressing multiple scales of sedimentary cyclicity previously interpreted to reflect annual to Milankovitch-scale forcing. Here we utilize X-ray fluorescence (XRF) core scanning and micro X-ray fluorescence (micro-XRF) scanning in combination with radioisotopic age data to evaluate a rock core record of laminated oil shale and carbonate mudstone from Utah's Uinta Basin, with the parallel objectives of elucidating the paleo-environmental significance of the sedimentary rhythms, testing a range of forcing hypotheses, and evaluating potential linkages between high- and low-frequency forcing. This new assessment reveals that the ∼ 100-μm-scale laminae—the most fundamental rhythm of the Green River Formation—are most strongly expressed by variations in abundance of iron and sulfur. We propose that these variations reflect changes in redox state, consistent with annual stratification of the lake. In contrast to previous studies, no support was found for ENSO or sunspot cycles. However, millimeter- to centimeter-scale rhythms—temporally constrained to the decadal to centennial scale—are strongly expressed as alternations in the abundance of silicate- versus carbonate-associated elements (e.g., Al and Si vs. Ca), suggesting changes in precipitation and sediment delivery to the paleo-lake. Variations also occur at the meter scale, defining an approximate 4 m cycle interpreted to reflect precession. We also identify punctuated intervals, associated principally with one phase of the proposed precession cycle, where Si disconnects from the silicate input. We propose an alternative authigenic or biogenic Si source for these intervals, which reflects periods of enhanced productivity. This result reveals how long-term astronomical forcings can influence short-term processes, yielding insight into decadal- to millennial-scale terrestrial climate change in the Eocene greenhouse earth. 

The Early Eocene Climatic Optimum (~53-50 Ma) represents the most recent episode of sustained greenhouse climate, during which the deep oceans were as much as 12°C warmer than today. The lacustrine Wilkins Peak Member of the Green River Formation (Wyoming, USA) is one of the premier locations to study this period of global warmth due to its rich terrestrial archive of climate dynamics, biology, and geomorphology. Using radioisotopic geochronology, cyclostratigraphy, sedimentology, and geochemistry, previous studies have leveraged this extensive record to evaluate the ancient lake system’s temporal evolution and response to climate. Much prior work on Green River Formation cyclostratigraphy, including that of Alfred G. Fischer, has focused on the evaluation of oil yield, a measure of organic richness. In this study, X-Ray fluorescence (XRF) core scanning of a basin center core, Solvay S-34-1, is used to produce a high resolution (5mm), continuous, multi-proxy elemental record of the complete Wilkins Peak Member, spanning 240 meters. This new geochemical assessment is a component of a larger multidisciplinary investigation that that is underway, including new magnetostratigraphic and radioisotopic geochronology. Elemental abundances for a range of measured elements, such as Si, S, Cl, K, Ca, Ti, Fe, Zn, Br, and Rb, are interpreted in terms of evaporitic, siliciclastic, and redox-sensitive sedimentation, and show variable responses at specific Milankovitch (eccentricity, obliquity, precession) and sub-Milankovitch time scales. Using this long high-resolution geochemical dataset of the Early Eocene Climatic Optimum, we consider potential linkages between Milankovitch forcing and sub-Milankovitch forcing, and plausible non-linear transfer functions that translate the astronomical insolation signal into the stratigraphic archive.